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Abstract 

We study the dependence on the temperature T of Casimir effects for a 
range of systems, and, in particular, for a pair of ideal parallel conduc- 
tors, h by I2, separated by a vacuum, a distance I3 apart, with h I3 
and h h- We study the Helmholtz free energy A T , combining Mat- 
subara's formalism, in which the temperature T appears as a periodic 
Euclidean fourth dimension of circumference It = hc/ksT with the semi- 
classical periodic orbital approximation of Gutzwiller. The latter was 
shown to be exact for parallel plates at T = 0. By inspecting the known 
results for the Casimir energy at T — in two cases of a rectangular 
parallelepiped, (h ^> I3 and h 3> h, and l\ 2> h and h 3> h), one 
is led to guess at the expression for A T of two ideal parallel conduc- 
tors, without performing any calculation. The result is a new form for 
A T , namely, A T = -(2Mi( 2 I 3 /ir 2 ) J27 3 = i J2n T =-oo L_4 ( n 3. n r), where 
L(ji3,riT) = [(2ri3i3) 2 + (jitIt) 2 ] 1 ^ 2 is the length of a classical periodic 
path on a two-dimensional cylinder section. This expression for the free 
energy is equivalent to others that have been obtained in the literature. At 
T — the semiclassical approach provides a finite and systematic approx- 
imation scheme in terms of classical paths that is useful when the normal 
modes of the cavity cannot be determined either explicitly or implicitly. 
Slightly extending the domain of applicability of Gutzwiller's semiclassi- 
cal periodic orbit approach, we here evaluate the free energy at T > 
in terms of periodic classical paths in a four-dimensional cavity that is 
the tensor product of the original cavity and a circle. The validity of this 
approach is at present restricted to particular systems. We also discuss 
the origin of the classical (h- independent) form of A T for T ~ 00. 

Key Words: Casimir energies, semiclassical theory, periodic orbits, finite tempera- 
ture 
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1 INTRODUCTION 



The best known of the Casimir effects is the force per unit area between two 
parallel ideal plates, separated by vacuum, a distance I3 apart, at a tempera- 
ture T = 00. Some fifty years later, that and other Casimir effects continue 
to intrigue both theoreticians and experimentalists. Indeed, interest has been 
growing. We list a few general references^, ||, |[ |f|. Of particular interest 
to us here are recent experimental reports by Lamoreaux^) and by Mohideen 
and colleagues [Q. Very recently the possible role of this effect in microelec- 
tromechanical systems has been examined^]. To avoid alignment problems one 
studies the force between a sphere of radius R and a wall, where the point on 
the sphere closest to the wall is at a distance I from the wall, with I <<C R. (The 
analysis of the force in this case is intimately connected to the force between 
walls |^, ^| and the experiment is simpler to perform. An experiment which 
might be even simpler to perform, because the force would be greater, is the 
force between a sphere and a segment of a spherical shell. We are here deal- 
ing with the focusing of virtual photons BErt. Focusing of virtual photons by a 
parabolic mirror has also been considered For the values of R and I stud- 
ied, the accuracy was such that a number of corrections to the case of ideal 
conductors must be considered. In particular, finite temperature effects must 
be accounted for. 

The corrections to the Casimir force for T > were first studied in a sem- 
inal paper by Lifshitz|12 in 1956. The paper studied dielectric slabs, and the 
temperature effects were a matter of dispute for years. The subject was first 
clarified in 1978 It was shown that Lifshitz's basic results for the tempera- 
ture corrections were correct, but that Lifshitz had erred in taking the limit as 
the dielectric constant goes to infinity when going from a dielectric to an (ideal) 
conductor. 

Our primary interest will be in rederiving the temperature effect for ideal 
walls but in a new form. We will utilize two sets of results, one by Matsubara||l4| 
for accounting for temperature effects in a general setting, and the second a 
formulation of Casimir effects at T — in terms of semiclassical extremal paths 
by the present authors^). The reformulation has the merit that the dominant 
paths are often rather obvious, with contributions that can be evaluated even 
when the proper frequency modes cannot be determined. 

We take this opportunity to mention two significant papers |fl6|| by Balian and 
Duplantier not referred to in fl5|| . The methods used in ||l6|| were quite different 
from those used in our paper. We had not fully appreciated the two papers at 
the time was published. 



2 IDEAL PARALLEL PLATES: A HEURISTIC 
APPROACH 

Matsubara showed pi that in studying a system in three-dimensional space in 
equilibrium at temperature T, it can be convenient to consider the system in a 
four-dimensional Euclidean space whose fourth dimension at finite temperature 
T > is compactified to a circle of circumference 

l T = hc/(k B T) . (1) 
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Thus, consider two ideal parallel conductors at temperature T, with an area 
A = separated a distance I3 in vacuum with I3 <5C h and I3 «. li. The force 
between the plates in this case is proportional to l\ and I2 and the force per unit 
area does not depend on either dimension. At finite temperature, the essential 
geometry in this case is that of a two-dimensional cylindrical surface of length 
I3 and circumference It- (See App. A) 

In our study Jll| of the plates at T = the analysis involved the determi- 
nation of the difference, Ap(E), of the spectral densities in the presence and 
absence of the conductors. Gutzwiller [^7| showed that this quantity is given in 
semi-classical approximation by the periodic paths of extremal classical action, 
in our case the periodic paths of extremal length. (Note that none of these paths 
has zero length and that a stationary periodic path must involve a minimum 
of two reflections.) In the case of parallel plates these are classical paths going 
back and forth between, and perpendicular to, the walls. 

The analogous analysis at finite temperature would involve studying the 
semiclassical Green function with a periodic dependence on the Euclidean time. 
The initial and final space-time points are still the same, but the Euclidean 
geometry is that of a cylinder. We prefer to proceed differently. 

For T > (and fixed volume), it is natural to work with the Helmholtz 
free energy A T rather than with £ T , the (Casimir) energy itself, and there is 
a natural choice of the form of A T . Thus, abandoning any pretense of rigor, 
we guess at the form of A T - there remains the determination of an over-all 
numerical factor - and justify our choice of A T afterwards. 

The semi-classical path approach has advantages that go well beyond pro- 
viding a simple physical picture of the origin of the Casimir effect at T > 0; 
as noted earlier, one can often determine the significant paths even when the 
natural frequency modes cannot be determined (neither explicitly nor, as in the 
case of the generalized argument theorem [p~8|, p"9[ , implicitly). 

We have 

A T — A T (h, c, l\, I2, h, It)- (2) 

At T — 0, we were concerned with the energy density. At T > 0, we will be 
concerned with the free energy density. In both cases we explicitly extract the 
volume I3A of the cavity. Further, on dimensional grounds, A T is proportional 
to he. We can therefore write 

A T = hcl 3 Ag(h,l T ) , (3) 

where g is proportional to the inverse fourth power of a length. Now, in |l5f it was 
shown that for a cavity with I3 <Kli,h, Sens is proportional to X^n 3 L~ (na), 
where the length L(ns) = 2713Z3 and 713 = ±1,±2,.... Further, for a paral- 
lelepiped with I3 « l\ and h « h, £c&s was proportional to J2n 2 Sn 3 ^ _4 ( n 2, ^3), 
where L(n2, nz) = [(2n2/2) 2 + {^n^l?,) 2 ] 1 / 2 and the sum is over all pairs of inte- 
gers (712,713), positive, negative, and zero, other than (0,0). £(713) is the length 
of a classical periodic path perpendicular to the walls and reflected 1 7T,3 1 times 
from each of them. Positive values of 773 refer to paths which start to the right, 
while negative values of 713 refer to paths which start to the left. L(n2, 773) is the 
length of a classical periodic path with | | reflections off each of the walls with 
dimensions l\ by I3 and 1 7x3 1 reflections off each of the walls with dimensions l\ 
by 12- For the case with Z3 « l\ and I3 « I2, T > 0, to which we henceforth 
largely restrict our attention, the knowledge that the space is enlarged by one 
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extra Euclidean dimension strongly suggests that A T can be obtained by simply 
replacing L(n2,n 3 ) by 

L{n 3 , n T ) = [{2n 3 l 3 ) 2 + (n T l T ) 2 ] 1/2 (4) 

This length is that of a classical periodic extremal path on the surface of a 
cylinder of height l 3 and circumference It- The path starts in a given direction 
and reflects |7T-3 1 times from each of the ends of the cylinder section (the ideal 
walls) and circles the cylinder in a particular direction \tit\ times. (The absence 
of a factor of 2 in the jitIt term of Eq. (f|) is due to the fact that we are 
considering extremal paths on a section of a cylindrical surface. See App. A) 
Negative values of n 3 refer to paths that start in the opposite direction, while 
negative values of tit refer to paths that circle the cylinder in the opposite 
direction. The contribution from the n 3 — path can be ignored since it gives 
an (extensive) contribution to the free energy that is proportional to the volume; 
it does not contribute to the (net) force on a plate. [See paragraph in Sec. 4 



above the paragraph which contains Eq. (21).] 
We thus guess that 

A T = KhcAl 3 L- 4 (n 3 ,n T ) , (5) 

(n3,UT) 

where the sum extends over all pairs of integers (n 3 , tit) with n 3 ^ 0, the length 
L(n 3 ,riT) is given by Eq. (^), and K is a numerical coefficient. We determine 
the value of K by demanding that A T reduce to 

tt 2 Tic A 

fcas -"720lf (6) 

for T ~ 0. But T ~ implies that It ~ oo which in turn implies that only 
tit — contributes to the sum in Eq. (JsJ) . For T~ we thus have 

A™ = KhcAl 3 2 f , * = ^ Khc f . (7) 

(The factor of 2 reflects the fact that we are now summing over positive values 
only of n 3 .) Comparison with Eq. (p|) gives K = —1/tt 2 and we obtain 

^ = ~^E £ L-\n 3l n T ). (8) 

713 — 1 tlt — — oo 

In App. B we show that Eq. (^) is indeed correct by comparing to a previous 
result. A proof from first principles, for a scalar field, is given in Sec. 5. The 
force per unit area at a temperature T is then given by 

^ = ii.T = _^V V 3(^ 3 ) 2 ~ (nrh) 2 (] 
A Adl 3 ^ 2^ n ^ [(n T l T ) 2 + (2n 3 l 3 YY ( > 

The somewhat complicated form of F T as opposed to the relatively simple form 
of A T shows that the free energy is not just a more fundamental object of study 
but also a more transparent one. Given A T , F T follows as in Eq. (^|). 



4 



3 THE HIGH-T AND LOW-T EXPANSIONS 

The free energy A T is readily obtained from Eq. (@) for z large or small, where 

z = 2l 3 /l T = 2l 3 k B T/(hc) . (10) 

Isolating the tit = contribution, the free energy of Eq. (g|) may be rewritten 
in the form 

A T ih)_ A (k B T)\ 

A ~ 720/1 tt 2 M 3 [ ' ■ 1 j 

From Eq. M) and the definition Eq. (0) we have that the dimensionless function 
A(z) is 

OO OO 

A (*) = E E ^ z + ™ 2 /^r 2 = a(i/z). (12) 

n— 1 m— 1 

[To check Eq. ([l2|), let z — > 1/z, and interchange m and n.] The reflection 
property |2(|, A(z) = A(l/z), is quite powerful. It relates the high- and low- 
temperature regimes, or, equivalently, the large and small l 3 regimes. The 
expansion of A(z) for z 1 is found by using the identity |2lf| 



— 2nna 



E ^2^2 = -W> = ^(wrcothC^) -l) = -_ + |L + ^e 

n—1 n=l 

- (13) 

Differentiating Eq. ( p~3| ) with respect to a 2 and then replacing a by toz one 
obtains, after multiplication by z 2 and summation over m, 



A(-) = E 



m= 1 



1 



2m 4 z 2 4to 3 z 2m 



— =— > (1 + 27rmnz)e~ 
rc 3 z ^-T 



ra=l 



f C(3) + f fV 2 ™* (1 + 2™*) ^m^ 3 , (14) 



180z 2 4z ,w 2z 

where C(3) = E™=i m ~ 3 = L2 02 . . .. Using the expansion of Eq. in Eq. @ 
gives for the free energy per unit area, for 2l 3 k B T > he, or 2/ 3 > It, 

A T ((3)k B T (k B T) M , Z T , _ W , T , 



.4 8vr/? Z T Z 



(1 + J^^-W^+of^e-s^A _ (15) 
3 i T i-i 47ri 3 \ 1 3 / 



Using the expansion of Eq. ([14]) with z replaced by 1/z in Eq. ( |TT| ) gives, for 
2l z k B T < he, or 2/3 < Z T , 

A T TT 2 hc kaT w . !T 2 knTU kaT , lo . ,111 „ {kaT /, 

-C(3)+ ^3— --yV (1+- ^-)e" 7rZr/i3 +0 -Z— e - 27TlT ' h 



A 720Zg 2tt/ 2 , sw 45/3 Z T Z 3 v ^ ; T Z 

(16) 

These low-T ( Eq. (jl6j)) and high-T ( Eq. (Jl5|) ) expansions of A T /A, first 
obtained by Sauer, by Mehra, and by Levin and Rytov[^2|, have also been 
reproduced in|l3[]. (The present derivation is somewhat simpler, to some extent 
because some of the earlier derivations were in a broader context, with dielectric 
walls or rather arbitrary surfaces. Note that Mehra|22 uses a convention for 
which an attractive force between the plates is taken to be positive.) See|l 
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for a physical explanation of the contribution proportional to X" 4 in the low-T 
expansion Eq. ( |l6| ) of the free energy. 

The expansions Eq. ([l6]) and Eq. ([l5]) are both exponentially good and the 
two limits are "dual" in the sense that since A(z) = A(l/z), they both follow 
from Eq. (|l4|). The relative error of these approximations is largest at the 
transition point, 2l%kBT = he, or z = 2l 3 /l T = 1. (The l 3 dependence of A T / A 
is not well defined at the transition point, since, with n arbitrary, a factor of 
(21 3 /It) u ~ 1 does not change the function appreciably near this point.) The 
exponential corrections all have the same sign and are no greater than ~ 2.3% 
at any value of l^ksT. Perhaps most significant experimentally is that the 
high-temperature limit is valid for I3 rather large compared to he/ (2ksT). This 
implies that the fall-off of the Casimir force at room temperature is proportional 
to Zg~ 3 rather than l^ 4 for separations ^3 rather large compared to 4 microns. 
Since the largest separation between the plates of early experiments was a few 
microns, it appears unlikely that a Casimir force proportional to l/l 4 could 
have been observed over the entire range. It should, however, be recalled that 
we have here considered an ideal metal rather than a real one. On the other 
hand, corrections due to finite conductivity diminish with increasing separation 
and are expected to be negligible compared to the temperature correction for 
separations exceeding a few microns || |23| . 



4 THE CLASSICAL FORM OF A T FOR T -> oo 

The fact that for parallel plates A T at sufficiently large T assumes a classical 
(^.-independent) form, as seen in Eq. JlE|), warrants some comment. We begin 
with the remark that A T in the high-T limit is not generally ^-independent. 
Thus, for a thin smooth closed surface, one finds Jig] 

A T = ~ak B T\n(k B Tl/hc) , (17) 

where the dimensionless quantity a is a function of ratios of the radii of cur- 
vature and I is an appropriately defined overall length scale. The logarithmic 
dependence of Eq. (|l7j ) on the overall length scale and ksT/hc implies that the 
change of the free-energy does not depend on he when all lengths are rescaled 
by the same amount, i.e., the "breathing mode" is always classical at high tem- 
peratures. 

The fact that the free energy is a pure power law for parallel plates makes 
that case a bit special. If we assume that the high-T limit of the free en- 
ergy of two parallel plates is classical, the power law dependence of A T fol- 
lows on dimensional grounds. (It follows from A T — A T (kBT, l\, I2, h) that 
A T = ksTli^gih), where g(l 3 ) must be proportional to I/Z3.) One may ar- 
gue for the existence of a classical high-T limit (and the absence of logarithmic 
corrections) from the contribution to the free energy of a single cavity mode of 
frequency u>k- Writing 



k k 
Aj for frequencies uj^ <k ksT/h becomes 



^+/r 1 ln(l-e-^ fc ) 



(18) 



A T k ~°° = ^+r 1 ln[/3? lWfc (l-^+0((/3^ fe ) 2 ))] ^fc B Tln(|^) , (19) 
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ignoring terms of order {3(hujk) 2 . Inserting the high-T expression Eq. ( fl9| ) 
in Eq. ([l8]) formally gives an asymptotic expansion of A T . But since the fre- 
quencies of a cavity generally can be arbitrarily high and the neglected terms 
are not small for all frequencies at a fixed value of T, this expansion is valid 
only in the presence of a cutoff. However, for the change of the free energies 
of the cavity in two different configurations the mode sum is finite if the two 
configurations can be connected adiabatically, as, for instance, by moving an in- 
ternal boundary arbitrarily slowly. An asymptotic expansion for large T of this 
difference does make sense. Under an adiabatic change of boundary conditions, 
the eigenmodes of the cavity change from lu^ to oj' k , but the quantum numbers of 
a mode are conserved. Using Eq. (^) one obtains for the asymptotic difference 
of the free energies when T <~ oo 



AA T = A' T — A 



k R T 



k 



(20) 



The sum in Eq. (g0|) converges if the difference of the free energies of the two 
configurations has a finite limit at T = 0, since a finite temperature does not 
induce new ultraviolet divergences. The expression Eq. (|2^) for the difference 
of the free energies at sufficiently high temperature (or better, for ksTl 3> he, 
where I is a relevant length scale of the problem) is classical and notably does 
not depend on h. In this limit the entropy S of the cavity is S — —A T /T 
and Eq. (|2(j) states that the difference in the entropies of a cavity at high 
temperatures approaches a constant that does not depend on the quantum scale 
h. (That the entropy of a bosonic gas in a cavity has a finite limit as T — ► oo 
is compatible with dimensional reduction, that is, the high temperature limit is 
described by Euclidean quantum field theory in 3 dimensions with h replaced 
byi/3@.) 

Although plausible, it is not possible to prove the existence of a finite high 
temperature limit using only classical arguments, because one can only compare 
the energy fcgT with the frequencies of the cavity by introducing a quantity with 
the dimension of Ti. If the asymptotic limit exists, and we have just argued that it 
does, it is independent of the value of h one has chosen to make this comparison. 
On dimensional grounds, the classical contribution to the change in entropy 
from a single (independent) mode can then depend only on the frequency ratio. 
Further, since the contribution to the (change in) entropy of (independent) 
modes is additive, classically the (change) in entropy due to a single mode is 
proportional to ln(w'/u>). Eq. ( pp| ) then follows as the classical change in the 
free energy of a cavity from the definition of ks- 

In studying the force between l\ by I2 plates at a separation ^3 one can 
evaluate the energy in the volume between the plates and in the infinite volume 
outside the plates and then determine the force by differentiation with respect 
to I3. It is a standard practice to avoid infinite volumes by considering a box 
of fixed volume]] l\ x I2 X L with an additional (movable) plate parallel to 
the l\ by I2 walls at a distance I3 from one of the walls [pi|. The total free 
energy of the box is the sum of the free energies of each of the two subvolumcs: 
Af ot (li, l 2 , L; l 3 ) — A T (li,l 2 , h) + A T (li,l2, L — I3). The force on the additional 
wall is —dAj ot /dl 3 . Extensive contributions to the free energy of the subsystems 
(proportional to their volume) do not contribute to the dependence of Aj ot on 

1 There should be no confusion with our usage of L to also denote path lengths. 
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l 3 and therefore do not result in a (net) force on the middle wall. Note that 
generally only the difference of the free energy of the box compared to the free 
energy of the box in some "standard" configuration, (say with the additional 
wall at ; 3 = L/2) i.e., Aj iff {h,l 2 , L; l 3 ) = Af ot {h,l 2 , L; h) - Af ot (l u l 2 ,L; L/2), 
is finite. 

With these arguments, we can, purely classically, obtain the high temper- 
ature limit of the change in the free energy, AA T , of a scalar field in a box 
of dimension l\ x l 2 x L, as we move an additional wall from a position at a 
distance l 3 from one of the l\ by l 2 sides to the middle of the box at L/2. Thus, 
from Eq. (|2p|), we have 

d 2 k ± 1 
(2tt) 2 2 ^ 

v ' Tt = — r 



AAl as =k B TA I V- 



n— — 00 



' u(7m/l 3 ,k±) uj(Trn/[L - l 3 ],k±) 
In — ; — - + In ; — - — 

Lj[2im/L,k_i_) ui(2irn/L,kj_) 

(21) 

where (jj(k 3 ,k±) = c-Jk 3 + k' 2 ,, with k 3 — irn/l 3 where l 3 is the appropriate 
distance from the wall. Using the integral representation 

Ma/6] = ijf~^[e-^ -e-^], (22) 



for the logarithms in Eq. (|2lj), the integration over kj_ can readily be performed 
on using polar coordinates. In addition we use the reflection property, 

00 00 

J2 e~ n2 * 2x = (™r 1/2 E e ~" 2/ * ' ( 23 ) 

n— — 00 n=— 00 

to rewrite the resulting sums over n. [ Eq. (|2^) and Eq. ( |l2| ) both relate the 
values of a function at the argument x and at the "reflected" argument 1/x. 
Eq. (p3|) is a particular case of Poisson's summation formula.] We then obtain 



v n— — 00 u 

(24) 

The n = contribution vanishes and X)^L-oo reduces to 2 J^Li. The integrals 
with n / converge term by term; with 



V- = ^. (25, 



A 5 / 2 2a 3 / 2 



we thereby arrive at 



AA Sass = -^^ I i- C (3)[i + - *,] , (26) 

which for L l 3 is the leading term of Eq. ( |l5| ) up to a factor of 2; the factor 
appears because we considered the scalar case here. [We strongly suspect that 
there may be a simpler classical derivation of Eq. (|2^) than the one given here.] 
A related problem is that of the Casimir energy, at T = 0, of two spinless 
charged particles in free space. There is a classical contribution in that case 
too, but it is not the leading term and the question arises if the term can be 
evaluated using classical theory. The answer in this case is also positive ]2l|. 
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An interesting point is that on dimensional grounds the absence of Ti leads 
to the absence of c, though we are studying electromagnetic waves. There is 
no contradiction; the classical Maxwell equations have meaning in the non- 
relativistic (c — > oo) limit|2(|. Thus, for example, in the dipole approximation 
the classical Maxwell equations give a c-independent radiative transition rate. 



5 THE SCALAR FIELD 

In App. B we prove the validity of Eq. (||) by making a connection with the 
results of|l3|. We here give an ab initio proof of that equation for a massless 
scalar field. The relation to the electromagnetic case is discussed in App. C (in 
a less broad context than that given, for example, in ref . [ll| ) . 

Thus, to substantiate some of the arguments of the main text, we consider 
a free massless scalar field <f> that satisfies Dirichlet boundary conditions 

<t>{ii = 0) = <j>{xi = k) = , (27) 

at the boundaries of a box with spatial dimensions l\ x I2 x . . . x Ip. The box 
serves as an infrared cutoff and guaranties a discrete spectrum. The dimension 
D of the box is arbitrary and can be used to dimensionally regularize the ultravi- 
olet behavior of subsequent expressions by analytic continuation to non-integer 
dimensions. Additional boundary conditions may be imposed on the field in- 
side the box. One could, for instance, demand that the field also vanish on a 
spherical shell located within the box, or, more appropriately for the Casimir 
force between plates, that it also vanish on a plane within the box. 

The free energy A T is related to the partition function Z by A T = — In Z, 
where (3 — l/(fcsT). Formally, Z is 

Z = Tre-P H , (28) 

where H = X)k^ k ( a ik ak + V^) is the free hamiltonian. ctk (a^) destroys (cre- 
ates) a quantum with quantum numbers k in the (possibly rather complicated) 
cavity. Crucial for the following is that H [ignoring any (weak) self-coupling 
of the scalar field] commutes with the number operators nk = a^ak no matter 
how complicated the boundary conditions on the field may be. H can therefore 
be diagonalized in a basis that enumerates the number of quanta in each cavity 
mode and has the eigenvalues E^n^ + 1/2) with each a positive integer 
or 0. Evaluating the trace of Eq. (p3) in this number basis, one has 

00 

In Z = ln(Q e- l3{nk+1/2)E ^} =-^2 ln[2 sinh(/3E k /2)] . (29) 

k nk=0 k 

Instead of computing the free energy or partition function directly, it will be 
convenient to first consider the energy at temperature T, 

Differentiating Eq. (PS), one finds 

£ T = £i?4+£e-^]. (31) 

k rri — 1 
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For £k > one can use Cauchy's theorem with a contour that runs just above 
the real axis and is closed by a large semicircle in the upper complex plane and 
write 

m\ + E e_m " £k ] - - lim + 1 / wi^^ 2 + E ^ > ( 32 ) 

2 ' £—0+ 7T J-oo+ie Et + A 2 ^— 

m=l ou-|-« K m=l 



up to terms that do not depend on E^ nor on j3\ in Eq. ( J32| ) we have dropped 
the contribution from the semicircular part of the contour. (For a sufficiently 
large radius O of the contour, this contribution is Q/tt up to terms that vanish 
for Q — ► oo.) The integral along the semicircular part of the contour thus 
does not depend on E^ nor T and can be absorbed in the overall normalization 
of the energy. Such constant contributions are irrelevant for the discussion of 
energy differences. With A = x + it and x real, the term on the right hand side 
of Eq. (|35), in square brackets, isp^] 



Q= [1/2+ e mP{x+lt) ] = l -cot{l3x/2)+n ^ 8{(3x - 2irn) + 0(e) . (33) 

m— 1 ~ n=— oo 



We can now evaluate the integral in Eq. (|32|). Apart from terms of order e that 
vanish in the limit e — > + , we can replace A by x. Because cot(f3x/2) is an 
odd function of x, the imaginary part of Q does not contribute to the integral. 
(This reflects the assumption that the cavity is ideal.) Formally, the energy of 
an ideal cavity at finite temperature is thus, from Eqs. (pL), (p2|), and (|33j), 



00 fO f) 00 

E E ^" 2 m = -|E E 

k n— — oo k 71= — oo 



n n e 



k n— — oo 



(34) 



where 



£(k, n) = [^ + (2t™//?) 2 ] 1/2 • (35) 

Since £ T = —(d/d/3)lnZ, Eq. j34| ) determines the partition function up to 
an overall normalization, or equivalently, the entropy up to an overall constant 
(which can be determined by demanding that the entropy vanish at T = 0). 
Eq. ( |34] ) formally shows that the partition function is proportional to the square 
root of the inverse determinant of the differential operator iff = h 2 — (hc) 2 -^-, 

where xt is the additional Euclidean coordinate in which the fields are periodic 
(with periodic length It) and h is the differential operator (single particle hamil- 
tonian) with eigenvalues E^. The positive eigenvalues of H4, the E(\t, n) given 
in Eq. ([55|), are the allowed energies of a massless scalar in the /owr-dimensional 
cavity constructed by augmenting every point of the original three-dimensional 
cavity by a circle of circumference It, that is, the 4-dimcnsional space is the 
tensor product of the cavity and a circle whose circumference is proportional 
to the inverse temperature. From Eq. ( pi[ ) and Eq. (^0|) we see that the free 
energy of the cavity can be written formally as 



^Lit y = r 1 E E ln£(k,n) = -/ dE'pl vity (E')HE') 



(36) 
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ith 

OO 

pLit y (£')=E E m*,n)-E'). (37) 



k n— — oo 



Eq. (37) shows that Pc av it y (E') * s the spectral density of the /owr-dimensional 



cavity constructed from the three-dimensional one in the above manner. Changes 
of p^ avity (E') due to (adiabatic) deformation of the cavity, including a change 
of temperature, can then be approximated semi-classically using periodic clas- 
sical rays of the /owr-dimensional cavity. This leads to a description of the 
free energy in terms of periodic paths in the /cwr-dimcnsional cavity which we 
exploited when we guessed at the form of the free energy in Sec. 2. 

To be somewhat more concrete, consider as an example a cubic cavity with 
sides of length I. The momentum components for the mode with energy E k are 
kih = riiTrh/l, and the energies in the corresponding four-dimensional cavity at 
finite temperature are 



E(k, n) = irhc 



h/2 



(38) 



as always, the length associated with the fourth (periodic) dimension, the analog 
of the lengths associated with the spatial dimensions, is It/2- 

Note that Eq. ([36]) is formal in two respects: the integral over the energy 
diverges due to the ultraviolet behavior of the spectral density, and Eq. ( p6| ) 
specifies the free energy only up to (a generally similarly divergent) contribu- 
tion proportional to T. Considerably more rigorous and physically relevant 
expressions can be obtained from Eq. ( |36| ) for changes in the free energy due to 
adiabatic deformations of the three-dimensional cavity. In physical terms the 
difference in the free energy is work that has to be done to change the boundary 
conditions for the field and should therefore be finite whenever it is physically 
feasible to do so. The ambiguity that Eq. (|3^) defines the entropy only up to a 
constant that does not depend on the temperature is removed by requiring that 
the entropy vanish at T = 0. 

Let us finally remark that for T — > and E^ fixed, 

oo 

lim Y /3- 1 \n[El + (27rn//3) 2 } 1 / 2 = E k /2, (39) 

/3— >oc ^ — 4 

n— — oo 

up to an infinite constant that does not depend on E k and therefore cancels 
in differences. [ Eq. (|39|) follows upon subtracting the -Ek-independent term 
ln[/i 2 + (27T7J//3) 2 ] 1 / 2 from the summand in Eq. ( |39| ) and noting that the so 
regularized sum defines a convergent integral in the limit [3 — > with the value 
(Ek — /i)/2. The dependence on the subtraction point \x can be absorbed in 
the overall energy normalization and cancels in energy differences.] At T = 
Eq. (^) thus gives the usual expression for the Casimir energy. The classical 
periodic orbits of the four-dimensional cavity in this limit are just the periodic 
orbits of the three-dimensional cavity, since any path that winds about the 
fourth dimension becomes arbitrarily long as T — > 0. Note that whereas the 
spectral density pj av i ty generally is proportional to the four-dimensional volume 
^rV'cavity the factor 1/It in Eq. ( ^6|) ensures that A T generally will only be 
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proportional to the three-dimensional volume of the cavity (which we extracted 
explicitly in Sec. 2). 

Although our starting point, Eq. ( |34| ) for the free energy, is the field the- 
oretic onc|L4|], the evaluation of this expression using the spectral density of a 
massless scalar in four spatial dimensions and, further, the semiclassical approx- 
imation to (changes in) this spectral density using classical periodic rays of the 
four-dimensional cavity is perhaps a bit unusual. This sketch of the relation 
between the free energy and the periodic classical rays of a corresponding four- 
dimensional cavity is the mathematical and perhaps even the physical basis for 
our guess of the free energy of two parallel plates in Sec. 2. 

6 DISCUSSION 

For ideal parallel plates, the semiclassical periodic orbital formalism gives a re- 
sult for £cas different in form but equivalent to that obtained by other methods. 
However, the approach provides a simpler physical picture, allows a simpler 
transition from £cas to A T , and gives a form simpler to evaluate numerically. 
It evaluates changes in the free energy due to changes in the temperature and 
in the cavity boundaries using the periodic classical paths in a four-dimensional 
cavity which is the tensor product of the original three-dimensional one and a 
circle. In a number of instances this amounts to replacing the contribution to the 
Casimir energy of each periodic classical path of length L 7 at T = by the sum 
of contributions of periodic classical paths of length L£(tit) = [L 2 + (titIt) 2 ] 1 ^ 2 
that, in addition, wind about the extra fourth dimension tit times. 

Furthermore, preliminary research indicates that this considerable simpli- 
fication may be present for a number of geometries, for scalar as well as for 
electromagnetic fields, and, at least for parallel plates, for dielectric walls. The 
simplification is a slight extension of the domain of applicability of Gutzwillcr's 
semi-classical periodic orbit approach, allowing the method to be used for a 
number of systems at finite temperature. Generalizing the result for parallel 
plates, it appears that the periodic classical paths of various lengths can be 
interpreted as dual variables to the eigcnfrcqucncics of a cavity: in rectangular 
cavities the two descriptions are related by Poisson's summation theorem. 
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APPENDIX A: FINITE TEMPERATURE AND 

THE INTRODUCTION OF A 
FOURTH DIMENSION 

As a vague indication of why the temperature can be treated by introducing 
a periodic fourth dimension, compare the standard operators exp(—iHt/h) in 
quantum theory and exp(— (3H) in statistical mechanics. One finds — it/%. 
The presence of i suggests periodicity, and we therefore drop the factor i and 
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replace f by a period r, giving (3 — t/K. For photons, we can replace t by 
It/c, leading to It — hc/ksT. For parallel plates, F/A must be proportional 
to Z1Z2, and the problem reduces to the determination of a function of l 3 and 
It- If one introduces periodic orbits between the plates, then the appropriate 
surface in l 3 and It is the surface of a cylinder of length l 3 and circumference 
It- Periodic classical paths involve lengths 2713/3, for 773 to and fro motions 
between the plates, and utIt, for circling 777- times around the cylinder. Since 
the coordinate system is Cartesian, the length of the shortest orbit that circles 
the cylinder 77^ times and goes back and forth between the plates 773 times is 
just Eq. (^). For a proper derivation of the introduction of a fourth dimension, 
see@ and Sec. 5. 



APPENDIX B: AGREEMENT WITH A PRE- 
VIOUS RESULT 

We now show that F T /A as given by Eq. (^) agrees with the force per unit area 
obtained previously @. Using Eqs.(3.8), (3.10), (3.11), (3.12), line 2 of (3.13), 
and (3.14) ofpf, we find 



El 
A 



where 



with 



c(a) = 
J [a) -- 



OO 

-3 c(nhc/(2k B Tl 3 )) 



n— — oo 



ATt 2 ah da 2 



J(a) 



sm(ay) 



d y- v 1 

ev-l 



(Bl) 



(B2) 



(B3) 



We obtain a different form for c(a) from that obtained in|l3| by expanding the 
denominator of Eq. (B2), 



00 „oo 

J(a) = Im / dye 
m=l J ° 



(ia—m)y . 



E 



a 



m z + a 2 



in agreement with Eq. ( |l3| ) . Performing the differentiation leads to 

/ s Sto 2 — a 2 

c («) = T^r 2^ 



2ir 2 l 3 ^ [to 2 + a 2 } 3 

rn—1 1 J 



(B4) 



(B5) 



Insertion of Eq. (|B5|) with a = nhc/(2k B Tl 3 ) = nl T /{2l 3 ) in Eq. (glj) repro- 
duces Eq. (0). 



APPENDIX C: ELECTROMAGNETIC VER- 
SUS SCALAR FIELDS 

We are primarily interested in electromagnetic fields and we therefore comment 
on the relationship between the free energy Aj c of a scalar field and the free 
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energy of an electromagnetic field, denoted, as elsewhere, by A T . The boundary 
conditions for scalar and electromagnetic fields are different, and further, for the 
latter case, both the condition div E = 0, which has no scalar field analog, and 
polarization must be accounted for; these three differences are not, of course, 
unrelated. Electrodynamics assumes different forms in spaces of a different 
number D of spatial dimensions and we will consider only D = 3. We will 
further limit our considerations to a rectangular box, with arbitrary dimensions 
h by h by h- 

As is well known, the condition div E = reduces the number of indepen- 
dent components of E to at most two and for a perfectly conducting box gives 
conditions on the quantum numbers m, n-i and 713 of the modes. There are two 
independent components if none of the rij vanish, one if one of the rij vanishes, 
and none if two of the rij vanish. In the periodic orbital approach one consid- 
ers classical ray's; in this case the condition div E = implies certain phase 
changes of each of the two independent polarizations upon reflection. 

We consider polarization effects for three cases, i) l\ ^> I3 and h ^ h, h) 
h 3> h and l\ 3> I3, and hi) arbitrary values of h,l2 and I3. For case i), the 
orbitals pass back and forth between, and perpendicular to, the two l\ by I2 
walls, the two states of polarization behave identically, and the effective electro- 
magnetic response function, 5g e i , is twice that of the relevant scalar response 
function Sg sc . For case ii), the situation is the same for orbitals which pass 
back and forth between the l\ by ^3 or the I2 by I3 walls. For orbitals which 
reflect off four walls, the pair of l\ by l 3 and the pair of I2 by l 3 walls, the 
polarization vector e is unchanged on reflection for e perpendicular to the plane 
of scattering; e is changed on reflection for e in the plane of scattering, but 
returns to its original value at the end of the (planar) periodic orbit, and once 
again we have, effectively, 5g e i = 2Sg sc . The situation is quite different for case 
iii), since periodic orbits that reflect off all six walls are not planar. [Periodic 
rays that reflect off all six walls cannot be planar, since no plane can cut all 
surfaces of the rectangular three dimensional volume.] To visualize one such 
non-planar orbit consider a classical ray that is close to and parallel to one of 
the spatial diagonals of the parallelepiped - not contained in any of the surfaces 
- and lies in one of the diagonal planes that contains the diagonal and two of 
the corners of the parallelepiped not connected by the diagonal. After reflection 
off all three walls at one of the corners, the ray, for symmetry reasons, is again 
parallel to the diagonal of the parallelepiped and in the same plane as before 
reflection at the corner. After reflecting off the other three walls at the diago- 
nally opposite corner, the path of the ray closes to a periodic one. If the initial 
ray is chosen very close to the diagonal, the length of this periodic orbit at the 
very least is close to twice the length of the diagonal of the parallelepiped, i.e., 
L(l, 1, 1) = [(2/i) 2 + (2l 2 ) 2 + {2h) 2 ] 1 ' 2 . The length of the periodic classical ray 
that reflects once off each pair of parallel walls in reality does not depend on how 
close to the diagonal the incident ray is chosen and thus is twice the length of 
the diagonal. The calculation would be much more complicated if that were not 
the case. More complicated non-planar periodic paths, that reflect a different 
number of times off each pair of parallel surfaces are obtained by changing the 
initial direction of the ray. For non-planar orbits there is a mixing of the two 
states of polarization and Sg e i generally is not just 25g sc ; the contribution to 5g e i 
of non-planar orbits is proportional to the trace of an orthogonal two-by-two 
polarization matrix. The existence of non-planar periodic orbits for case iii), 
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and the consequent non-diagonality of the polarization matrix, make less sur- 
prising the fact that £cas in case iii) depends p8| not only on the 3-dimensional 
path lengths, 

L(m,n 2 , na) = [(2Zim) 2 + (2l 2 n 2 ) 2 + (2Z 3 n 3 ) 2 ] 1/2 , (CI) 

but also on the one-dimensional path lengths L(rij) = 21^ for i = 1, . . . , 3. The 
presence of one-dimensional paths can also be associated with the restriction of 
the quantum numbers due to the boundary conditions. Thus for scalar fields 
satisfying Dirichlet boundary conditions in a rectangular cavity - the eigenfunc- 
tions are products of sines - none of the m quantum numbers vanish. Finally, 



we note that we stated in| 15 that the semiclassical approach was exact if any 
of the ratios formed from 1%, l 2 and I3 was arbitrarily large, as is the case for i) 
and ii) but not for iii). The restriction was necessary because we had assumed 
that Sg e i = 2Sg sc . 

We remark, parenthetically, that the four dimensional space characterized by 
li, l 2 , £3 and It, with three spatial dimensions and one periodic dimension (with 
a circumference proportional to the inverse temperature), is not equivalent to 
electrodynamics in four spatial dimensions; div E 7^ for the latter space. 
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